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We perform three-dimensional relativistic hydrodynamical calculations of neutron star mergers to 
assess the reliabihty of an approximate treatment of thermal effects in such simulations by combining 
an ideal-gas component with zero-temperature, microphysical equations of state. To this end we 
compare the results of simulations that make this approximation to the outcome of models with 
a consistent treatment of thermal effects in the equation of state. In particular we focus on the 
implications for observable consequences of merger events like the gravitational-wave signal. It is 
found that the characteristic gravitational-wave oscillation frequencies of the postmerger remnant 
differ by about 50 to 250 Hz (corresponding to frequency shifts of 2 to 8 per cent) depending on the 
equation of state and the choice of the characteristic index of the ideal-gas component. In addition, 
the delay time to black hole collapse of the merger remnant as well as the amount of matter remaining 
outside the black hole after its formation are sensitive to the description of thermal effects. 
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I. INTRODUCTION 

High-density matter above nuclear saturation density 
as it is present in the cores of compact stars usually oc- 
curs at temperatures that are low compared to the Fermi 
energy. Only during the formation process in a stellar 
core collapse accompanied by a supernova explosion are 
temperatures of several 10 MeV reached, where, however, 
because of neutrino losses the temperature drops below 
1 MeV within the first minute. Higher temperatures can 
only be reached again in neutron star binary systems at 
the end of a tens to hundreds of Myrs lasting inspiral 
phase, where due to gravitational-wave emission the or- 
bit shrinks and a merging of the two binary components 
takes place (see [l[ for a review). 

About ten binary neutron stars are known in our 
Galaxy, while the actual population is assumed to be 
much larger due to selection effects Although no 

direct observations of the unavoidable merging processes 
have been made, it is clear from simulations that the tem- 
perature rises up to some 10 to 100 MeV in the merged 
objects (see e.g. Q). It is also known that the addi- 
tional pressure support due to thermal effects influences 
the structure and development of the merger remnants 

The modeling of neutron star mergers is limited by 
the incomplete knowledge of the equation of state of 
high-density matter. A large variety of EoSs have been 
proposed based on various theoretical descriptions and 
matched to various observational data of nuclear matter 
(see e.g. Neutron stars described by these EoSs show 
a big spread of the stellar properties, e.g. of the compact- 
ness and the maximal mass that can be supported against 
gravitational collapse @. (The astronomical observa- 
tions of theses properties are not conclusive to constrain 
the EoSs significantly [l3|, but see Most EoSs con- 

sider the ground state of matter, i.e. at zero temperature 
and in equilibrium with respect to weak interactions (in 
the following referred to as "cold" EoSs). However, in 



the case of compact star coalescences thermal effects are 
of significant influence 0-01 • Currently, there are only 
two microphysical EoSs that include nonzero tempera- 
ture effects and that are in use for such kind of studies 
(see [l for a review) : the model by (l2j , which we will re- 
fer to as LS in the following, and the calculations by [l3j . 
which will be called Shen. In addition, an approximate 
treatment of thermal effects in combination with any cold 
microphysical EoSs has been employed e.g. in jl4| - |l7l |. in 
particular in current fully relativistic studies with mi- 
crophysical EoSs. However, the reliability of this ideal 
gas-like extension has not yet been explored. 

It is the purpose of this paper to investigate the va- 
lidity of this approximate description of thermal effects. 
An understanding of the limitations of this approach is 
important for the interpretation of observationally rele- 
vant results of simulations adopting this scheme. This in- 
cludes especially the gravitational- wave signals of merger 
events, which may be detectable with current and up- 
coming grav itational- wave detectors like LIGO [3 and 
VIRGO [1^. Such experiments have some potential to 
constrain the high-density EoS once measurements be- 
come available. These attempts will have to rely on the 
results of simulations with various theoretical prescrip- 
tions of the EoS using partially th e ap proximative tem- 
perature implementation (e.g. [14l4l7l |). which makes it 
important to estimate the influence of this simplified de- 
scription. Furthermore, neutron star mergers are specu- 
lated to be the progenitor systems of short gamma-ray 
bursts, and the modeling of the conditions for producing 
such gamma-ray bursts may also depend on the particu- 
lar treatment of thermal effects 0, [13] ■ 

First we briefiy review the relevant earlier findings, 
which our analysis is supposed to supplement. Our meth- 
ods are introduced in Sect. IIIII and the properties of the 
EoSs used in this study are discussed in Sect.|TVl Sect.lVl 
presents the temperature phenomenology in merger simu- 
lations. In Sect. lVII we investigate the influence of the ap- 
proximate treatment of thermal effects on observational 
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features of the merging process, and Sect. I VIII summa- 
rizes our conclusions. 



II. REVIEW OF PREVIOUS STUDIES 

Different approximations of microphysical EoSs have 
been compared in previous studies, which we briefly sum- 
marize here to set the stage for our analysis. We focus 
on those works that examine and demonstrate the impor- 
tance of including thermal effects in neutron star merger 
simulations. 

The simplest approach to describe neutron star matter 
is the use of the polytropic relation P = Kp^ , which deter- 
mines the pressure P as dependent on the rest-mass den- 
sity p with constant k and F. Employing this barotropic 
function during the evolution of a neutron star merger 
corresponds to the extreme case of perfectly efficient 
cooling, since this EoS does not allow for shock heat- 
ing ("iscntropic case"). This means that thermal energy 
is taken away from the matter instantaneously, and in 
fact the energy equation does not need to be evolved. As 
an alternative, one can use an "ideal-gas" ansatz, where 
the pressure is given by P = (F — l)ep with the internal 
specific energy e. In this case the energy is conserved also 
when shocks develop. This approach is particularly in- 
teresting if the initial data are constructed also by a poly- 
tropic relation and a comparison between iscntropic and 
ideal gas simulations is made. This allows one to judge 
the importance of thermal effects. Both EoS treatments 
were compared in where the lower pressure support 
in the iscntropic case led to an earlier gravitational col- 
lapse of the merger remnant compared to the ideal-gas 
calculation. Furthermore, the matter remaining outside 
the black hole formed a thin disk-like structure, while 
the ideal-fluid EoS resulted in a vertically inflated torus. 
This finding can also be explained by the additional pres- 
sure support provided by the nonzero temperature, which 
"inflates" the low-density material in the outskirts of the 
remnant. 

In Rcf. [H it was investigated to which extent an an- 
alytical ideal-gas EoS can be used to reproduce the fea- 
tures of a neutron star merger simulation with a micro- 
physical EoS. The initial stars were set up as polytropes, 
with At and F fixed in a way that the relation between 
central density and gravitational mass of isolated neutron 
stars matched the relation of the microphysical EoS in a 
wide range. It was found that the simplified treatment 
yields a remnant structure and a density profile simi- 
lar to the one obtained with the microphysical EoS. The 
gravitational- wave emission however, was only in quali- 
tative agreement plj . and not all relevant features were 
compatible between the simulations. Moreover, these 
conclusions were based on only one binary configuration, 
and it remains unanswered whether the similarities also 
hold for a larger variety of models. Furthermore, the 
study was performed for only one microphysical EoS and 
it should be explored whether this approximation works 



also for other microphysical EoSs. 

In [11, 0] it was analyzed whether the inclusion of ther- 
mal effects is at all important for the remnant properties 
and observable quantities in neutron star merger sim- 
ulations when microphysical EoSs are used. The com- 
parisons presented there are closely related to the afore- 
mentioned in 0. The data of a temperature-dependent 
microphysical EoS were implemented in a way to impose 
perfectly efficient cooling, i.e. only the zero-temperature 
sector of the EoS was taken into account. This model was 
compared with simulations employing the same micro- 
physical EoS, but using its full temperature dependence. 
As in the major differences are explained by the lower 
thermal pressure contribution in the zero-temperature 
calculations. While in the case of the temperature- 
dependent simulation an inflated, "hot" torus was ob- 
tained, the omission of thermal effects led to a disk with 
a small vertical extension and a higher density. In addi- 
tion, the "cooling losses" in the zero-temperature mod- 
els resulted in more compact merger remnants, which 
shifted the gravitational-wave frequencies in the post- 
merger stage to higher values. Moreover, it also increased 
the amount of matter remaining outside the black hole af- 
ter the gravitational collapse of the central object [7l[2l|. 

In summary, the studies mentioned here show that the 
consideration of nonzero temperature effects is generally 
important in simulations of neutron star mergers. There- 
fore, in the following we will analyze a more sophisticated 
approach to include thermal effects in an approximate 
manner in simulations with microphysical EoSs. It was 
originally introduced in [2^ and has been applied in the 
context of stellar core-coUapse (e.g. [l^]) and in neutron 
star merger studies (e.g. jI44l7| V The ansatz provides an 
additional pressure contribution Pth = (Eth ~ l)eth due 
to thermal effects, which is proportional to the thermal 
energy density eth assuming an ideal-gas like behavior 
with a constant "ideal- gas index" Fth- The value of Fth 
cannot be unambiguously determined from physics argu- 
ments but was chosen more or less ad hoc. It will be 
argued in Sect. IIVI (see Fig. [2|) that assuming a constant 
ideal-gas index the applicability of the ansatz is at least 
questionable and requires an investigation. Furthermore, 
wc will explain in Sect. IIIII that for barotropic EoSs this 
approximate description implies an unphysical variation 
of the electron number, which affects spuriously the EoS 
properties and thus may have an impact on the dynamics 
of a neutron star merger. 



III. METHODS AND SETUP 

We perform simulations of neutron star mergers using 
a general relativistic Smoothed Particle Hydrodynamics 
scheme, which employs the conformal flatness approxi- 
mation of general relativity [2^ . The description of the 
code can be found in [H, [25| . Various microphysical EoSs 
are implemented in our model to close the set of hydro- 
dynamical equations. The EoS provides the pressure P 
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and the specific internal energy density e as a function 
of the rest-mass density p, the temperature T and the 
electron fraction Y^. The rest-mass density, the energy 
density and the electron fraction are given by the hydro- 
dynamical evolution. We assume the advection of the 
initial electron fraction and neglect neutrino production, 
which leads to the simple evolution equation ~ 0. 
A first look-up in the EoS determines the temperature 
as a function of the energy density for the given p and 
Ye by a numerical inversion procedure. Then the pres- 
sure is found from the EoS as dependent on p, Yq and T. 
Both operations involve numerical interpolation schemes, 
because the EoS is given in the form of a table. 

For this study we employ the Shen EoS and the LS 
EoS. These EoSs provide the full temperature depen- 
dence and therefore, they can be used for a comparison 
to quantify the deviations between a consistent and an 
approximate inclusion of thermal effects. 

The approximative description of thermal effects can 
be applied for any cold, microphysical EoS. For this rea- 
son it was included in simulations of neutron star merg- 
ers. The zero-temperature EoSs provide an unambiguous 
relation between the rest-mass density and the pressure, 
as well as the rest-mass density and the energy density. 
Typically, these functional dependences are given in the 
form of a table. For barotropic EoSs that arc commonly 
used, the electron fraction is fixed by equilibrium with 
respect to weak interactions, i.e. the evolution of the ad- 
vection equation = used in consistent models is 
abandoned [ssj . 

In order to include temperature effects one splits up 
the pressure and the specific internal energy into a cold 
and a thermal part: 

P = Pcold + Pth, (1) 

e = Ccoid + eth- (2) 

The cold contributions are taken from the tabulated, mi- 
crophysical, zero-tempcrature EoS and are functions of 
the rest-mass density p only. The rest-mass density p 
and the specific internal energy e are given by the hy- 
drodynamical evolution. Then, for a given e and p one 
defines the thermal part of the specific internal energy by 

eth = e - ecoid(/o), (3) 

where the cold contribution is taken from the cold mi- 
crophysical EoS table determined by p. Assuming an 
idcal-gas like behavior the thermal contribution to the 
pressure is set to 

Pth - (Fth - l)peth, (4) 

with the ideal- gas index Fth, which is assumed to be con- 
stant for all p and e. Now one obtains the pressure from 
Eq. dl]), which closes the hydrodynamical equations. By 
means of this scheme arbitrary barotropic EoSs can be 
implemented in hydrodynamical simulations consistently 
with the first law of thermodynamics. 



Besides the approximate treatment of temperature ef- 
fects, a further, less obvious, approximation is implied 
by this implementation, when used with barotropic EoSs. 
The cold EoSs are constructed for matter in neutrino-less 
beta-equilibrium, i.e. for a neutrino chemical potential of 
zero. This means that the pressure and the energy den- 
sity arc functions of p only, and the explicit dependence 
of both on Yo is neglected. Therefore, when matter is 
compressed or expanded, Y^ changes unphysically due to 
the imposed condition of beta-equilibrium (because for 
a barotropic EoS Y^. in neutrino-less beta-equilibrium ef- 
fectively is a function of p only). However, in reality, 
matter that was initially in beta-equilibrium, would be 
out of equilibrium after compression or decompression, 
because the composition changes on a longer than the 
dynamical timescale [33 |. Hence, the instantaneously en- 
forced adjustment to beta-equilibrium leads to spurious 
effects in the dependent quantities of the EoS like the 
pressure and the energy density. 

For the sake of a comparison we construct cold EoSs 
from the tables of the LS and the Shen models, i.e. we 
keep only the information about the microphysical prop- 
erties at T = and beta-equilibrium. In fact, we use the 
barotropic relations that are also employed for computing 
the initial data of our fully consistent simulations. 

For both EoSs (LS and Shen) we perform three runs 
each: one using the full temperature dependence of the 
EoS, and two simulations where the approximate treat- 
ment of thermal effects is employed in combination with 
the cold EoS (see Tab. IJ). One calculation is done with 
Fth = 1.5, the other with Fth = 2. The choices for Fth 
wiU be discussed in Sect. IIVI 

For every model we simulate the merging of two non- 
spinning stars with gravitational masses in isolation of 
1.35 Af©. This setup can be considered as a canonical 
configuration that is predicted to be very abundant in 
the neutron star binary population [26| . Furthermore, for 
both EoSs this set of binary parameters results in the for- 
mation of a so-called hypermassive postmerger remnant, 
a differentially rotating object that collapses to a black 
hole after angular momentum redistribution on typical 
timescales of some 10 milliseconds. Higher total binary 
masses would cause the prompt formation of a black hole 
shortly after the stars have come into contact. Since 
temperature effects are important in particular during 
the formation and the evolution of a hypermassive post- 
merger remnant, we focus on the aforementioned equal- 
mass binaries with moderate neutron star masses. Our 
analysis is constrained to this limited, but representative, 
set of models because of the high demand of computa- 
tional resources. In particular, we found it better to com- 
pute until black hole formation if possible for the chosen 
configurations, rather than conducting more models for 
a shorter time. Therefore, our investigation should be 
considered as a matter-of-principle study rather than a 
detailed and full survey of uncertainties in the parameter 
space (masses, spins, EoS) of binary setups. 

Two additional simulations with different resolutions 
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TABLE I: Models discussed in the text. The osciUation fre- 
quency of the postmerger remnant is denoted by /peak. Tdeiay 
is the delay time between the merging and the formation of a 
black hole. Note that in the cases where the black hole col- 
lapse does not occur during the simulation, we give a lower 
limit for Tdeiay. This value is determined by the finite simula- 
tion time and may differ significantly from the true value (see 
[2?!'). The estimated mass of the torus remaining outside the 
black hole after the gravitational collapse is listed in the last 
column. For the simulations based on the Shen EoS we give 
the amount of matter fulfilling the torus criterion 16 ms after 
the merging, because for these models a black hole does not 
form within the simulation time. "LS 140K" and "LS 340K" 
refer to runs with lower numerical resolutions using the full 
LS EoS table. To estimate resolution effects these simulations 
should be compared to the entries in the first row. 
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EoS 


/peak [kHz] 


Tdeiay [ms] 


Mtorus [M0] 


LS full table 


3.24 


20.8 


0.064 


LS Pth = 1.5 


3.38 


10.3 


0.045 


LS Fth = 2 


2.99 


21.0 


0.054 


Shen full table 


2.21 


>18.1 


0.075 


Shen Tth = 1.5 


2.27 


>16.8 


0.084 


Shen Tth = 2 


2.16 


>18.0 


0.074 


LS 140K 


3.22 


17.5 


0.126 


LS 340K 


3.25 


18.4 


0.073 



are performed for the LS EoS using the consistent tem- 
perature description to assess the resolution dependence 
of the quantities we are interested in. For the character- 
istics of a merger discussed below, Tab. |T] gives the re- 
sults of runs with about 140,000 and 340,000 instead of 
our standard resolution of about 550,000 SPH particles. 
They show a strong dependence of the estimated torus 
mass on the resolution (deviations to the run with the 
standard resolution of about 97 per cent for the lowest 
resolution and of about 14 per cent for the intermediate 
resolution). The delay time to the black hole formation 
depends moderately on the resolution, while the oscilla- 
tion frequency of the postmerger remnant is less sensitive 
to the particle number. However, we stress that for each 
set of models that will be discussed in the following, based 
either on the Shen EoS or the LS EoS we use exactly the 
same initial data and the standard resolution. Therefore, 
deviations among these calculations are entirely caused 
by the different treatments of thermal effects. 



IV. EOS PROPERTIES 

The dependence of the pressure on the temperature 
is shown in Fig. [T] for both EoS models used in this 
study. Two constant values of the rest-mass density, 
p = 3 • lO^^g/cm^ (black curves) and p = 6 • lO^^g/cm"^ 
(red curves) , were chosen corresponding to about one and 
two times nuclear saturation density. For low tempera- 
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FIG. 1: Pressure as a function of the temperature for rest- 
mass densities of p = 3 ■ 10^''g/cm^ (black curves) and p = 6 . 
lO^^g/cm'^ (red curves). The pressure is shown for a constant 
value of the electron fraction as determined from the condition 
of beta-equilibrium at T = for the given densities. The solid 
lines correspond to the Shen EoS, while the dashed curves 
display the pressure for the LS EoS. 



tures (below 5 MeV) the pressure is nearly constant as a 
consequence of the degeneracy of nuclear matter at these 
densities. One also recognizes that the Shen EoS pro- 
vides a higher pressure at the same density. For the de- 
pendence of the pressure on the rest-mass density and 
an overview of the mass-radius relations of neutron stars 
described by the soft LS EoS and the stiff Shen EoS we 
refer the reader to Figs. 1 and 2 in Q. As a consequence 
of these EoS properties the stars described by the Shen 
model are in general less compact. In addition, pro- 
vides information about the influence of thermal effects 
on the stellar structure: for temperatures of some 10 MeV 
the stellar radii are a few kilometers bigger. 

For a given p, T and one can compute the local 
effective "ideal-gas index" Fth of the EoSs. According to 
Eqs. dl]), ^ and (g]) one finds 



rth(p,T,re) = — + 1 
peth 



P{p,T,Y,)~P{p,Q,Y,) 
p[e(p,T,yc)-e(p,0,re)] 



f 1. 

(5) 

This thus defined ideal-gas index Fth is shown as a func- 
tion of the rest-mass density in Fig. [2] for typical tem- 
peratures occurring during the merging of neutron stars 
(see Sect. |V|. The variation of the displayed values of 
the ideal-gas index reveals that the choice of a constant 
value of Fth can only be a very crude approximation and 
requires a validation, which we aim at with this study. 
Since we would like to assess the validity of this ansatz in 
neutron star merger simulations, we adopt the values for 
Fth previously used in the corresponding literature [l^ - 
\VA [22| . i.e. Fth = 1.5 and Fth = 2. One should keep 
in mind that the dynamics of the merger are determined 
mainly by the supranuclear EoS. Hence, the ideal-gas 
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FIG. 2: "Ideal-gas index" Fth as a function of the rest-mass 
density for the Shen EoS and the LS EoS for typical tem- 
peratures. See Eq. [S]for the definition of Fth- The electron 
fraction for a given density is determined from the condi- 
tion of neutrino-less beta-equilibrium at T = 0. 



index should primarily reproduce the thermal EoS be- 
havior in this density range. From Fig. [2] one might eon- 
elude that the ehoiee Fth = 2 is slightly too high, thus 
it could overestimate the thermal pressure contribution 
Pth, while Fth = 1-5 might yield a too low pressure sup- 
port. One also recognizes that the ideal- gas index is lower 
the higher the temperature is. In comparison to the LS 
EoS, the Shen EoS has higher values of Fth. 

Furthermore, for densities below ~ lO^^g/em'^ the 
ideal-gas index Fth approaches 4/3. This is a consequence 
of the fact that in this density regime the pressure is pro- 
vided mainly by an ideal gas of ultra-relativistic electrons 
and photons. 



V. SIMULATIONS 



Moreover, material is shed off from the merger remnant 
to form a dilute halo around the hypermassive object. 
The halo matter is shock heated by the interaction with 
the surface of the deformed, dense, rotating remnant, 
leading to temperatures of 20 to 30 MeV (upper right 
and lower left panels). Because of the differential rota- 
tion the shock-heated, hot material from the clash of the 
two stars is spread within the central object (lower left 
panel). After some milliseconds an approximately sta- 
tionary configuration is reached with a nearly axisym- 
metric temperature distribution (lower right panel). Re- 
markably, the outer part of the central core remains at 
a slightly lower temperature than the inner core and the 
surrounding medium (barely visible), because it is not as 
strongly compressed and it is not subject to shock heat- 
ing. During the whole evolution after the merging the 
temperature of the supranuclear matter forming the cen- 
tral object is above 20 MeV reaching even more than 
60 MeV at the center. From this analysis it is clear 
that the effects of nonzero temperature should be con- 
sidered as important for the dynamics and the shape of 
the postmerger structures (see Fig. [1] and the summary 
in Seet.ini). 

Basically the same picture arises for the simulation 
with the Shen EoS, where the merging proceeds similarly 
to the LS case. However, the temperatures are lower. 
When the remnant has settled to an axisymmetrie config- 
uration the central temperature is about 30 MeV, while 
the outer remnant parts have temperatures between 10 
and 20 MeV. This can be understood from the fact that 
the stiffer EoS leads to a less compact remnant structure, 
thus a less extreme compression of the matter. 

For the simulations using the approximate description 
of thermal effects, a relation between the thermal en- 
ergy and the temperature is not provided by the EoS 
and therefore temperature information cannot be given 
on this level. Only by introducing an additional assump- 
tion about the nature of the hot matter, for instance 
considering it as a gas of nucleons, one can estimate the 
temperature (see e.g. [HI). We refrain from doing this 
here. 



We evolve neutron star binaries with our hydrodynami- 
cal code through the last phase of the inspiral, the merg- 
ing and the postmerger stage. These different phases 
lead to distinctive temperature distributions in the or- 
bital plane of the merger with the LS EoS using the full 
temperature description (Fig. [3]). When the stars or- 
bit around each other, numerical dissipation heats the 
stellar interior to some MeV (upper left panel). As can 
be seen in Fig. [1] this temperature increase can be ex- 
pected to be dynamically unimportant. Then, during 
the merging, shock heating produces temperatures up to 
about 100 MeV at the contact layer between the stars, 
where also Kelvin-Helmholtz shear vortices develop (up- 
per right). The subsequent compression leads to an av- 
erage temperature of about 20 MeV in the dense core 
of the central object (upper right and lower left panels). 



VI. RESULTS 

In this section we discuss the consequences of an ap- 
proximative ideal-gas description of thermal EoS contri- 
butions for observable features of binary neutron star 
mergers. 

A. Gravitational waves 

The overall dynamics of the Fth-runs are similar to the 
corresponding simulations using the full temperature de- 
pendence. However, in order to judge the quality of the 
approximate treatment of thermal effects, one should in- 
vestigate the implications of the different approaches for 
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FIG. 3: Temperature in MeV in the orbital plane of a merging NS binary with 1.3 Mq and 1.35 Mq components for the LS 
EoS using the fully consistent temperature implementation. The rotation is counter-clockwise. Note that the temperature is 
given on a logarithmic scale in order to resolve all relevant features. The plots were created with the visualization tool SPLASH 



observable aspects. Therefore, we focus in our analysis 
on the gravitational- wave signal, which is the most direct 
probe of neutron star collision events. 

As can be seen in Fig. [3l temperature effects can 
only be dynamically important in the merging and post- 
merging stages, while they are irrelevant during the 
inspiral of the binary components. It is known that 
the gravitational-wave emission of the deformed, os- 
cillating, hypermassive remnant depends sensitively on 
EoS properties 0, [2l|. And it was also shown previ- 
ously that the inclusion of thermal effects has an im- 
pact on the gravitational- wave signal of the postmerger 
phase [1,0, ill. 

The emission of gravitational radiation from the post- 
merger remnant is dominated by a narrow frequency 
band around a frequency /peak, which is visible as a 
prominent peak structure in the spectrum (see Fig. ^ . 
Hence, this peak frequency /peak can be used as a quanti- 



tative feature of the gravitational- wave signal in order to 
characterize the emission during the postmerger phase. 

In Fig.|3]one can read off the peak frequencies of the LS 
models from the gravitational- wave luminosity spectrum, 
which is given by dEcw/df = 2Tr'^ D'^ f {\h+\^ + \h^\'^). 
Here D is the distance from the source, which cancels out 
when multiplied with the Fourier transformed waveforms 
h for both polarizations averaged over all emission direc- 
tions. The amplitudes of the gravitational- wave train are 
computed by means of a quadrupole formula that takes 
into account post- Newtonian effects 0, [2^ . 

With the LS EoS using the full temperature depen- 
dence, the peak frequency of the postmerger oscillation 
is located at 3.22 kHz (black curve in Fig.U]). In compar- 
ison, the peak frequencies of the simulations employing 
the approximate scheme are about 200 Hz lower for the 
Tth = 2 model and about 200 Hz higher for the Tth = 1.5 
run (sec Tab.|T|. To have a reference for judging the im- 
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FIG. 4: Direction and polarization averaged gravitational- 
wave luminosity spectra for symmetric binaries with Mi = 
M2 = 1.35 M© described by the LS EoS using the full tem- 
perature dependence (black curve) and an approximate de- 
scription of thermal effects with Fth = 1.5 (red curve) and 
Fth = 2 (green curve). Note that our simulations start only a 
few orbits before the merging of the stars. Hence, power from 
the preceding inspiral phase at lower frequencies (< 1 kHz) is 
missing in this plot. 
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FIG. 5: Enclosed gravitational mass within an ellipsoid 
with the semiaxes a = x, b = x and c = x/2 for the 
1.35 Mq-I-1.35 Mq binaries about 8 ms after the merging. 
Shown are results with the LS EoS using the full temperature 
dependence (black curve) and an approximate description of 
thermal effects with Fth = 1-5 (red curve) and Fth = 2 (green 
curve). Note that within general relativity the gravitational 
mass is only defined in isolation. Here Mgrav denotes the en- 
closed contribution to the ADM mass, neglecting the extrinsic 
curvature terms. Using the enclosed rest mass yields a similar 
result. The distance x is given in isotropic coordinates. 
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FIG. 6: Same as Fig. H but for the Shen EoS. 



portancc of such difFcrcnces, we note that for instance 
for a merger of a binary with 1.2 Mq and 1.35 Mq neu- 
tron stars described by the full temperature-dependent 
LS EoS the peak frequency is 3.04 kHz (see 0). This 
value differs also by about 200 Hz from our symmetric 
standard case model. 

The discrepancies between the models with the differ- 
ent treatments of thermal effects can be understood as 
a consequence of the different remnant structures. More 
precisely, the resulting compactness of the remnants is 
the crucial property determining /peak- The more com- 
pact the objects are, the higher frequencies are obtained. 
In this context we refer the reader also to 0, where the 
dependence on the compactness was discussed. Since the 
Fth = 1.5 implementation of thermal effects leads to a 
lower total pressure support, the remnant is more com- 
pact compared to the case with the full EoS. In contrast, 
Fth = 2 overestimates the thermal pressure contribution 
and consequently, the central object is less compact. 

For isolated stellar objects the compactness is defined 
as C = GM/ [c^R), with the gravitational mass M and 
circumferential radius R of the star. G and c are the 
gravitational constant and the speed of light. For neu- 
tron star merger remnants, however, this definition of the 
compactness is not useful, because the remnants do not 
have a well-defined surface (see Fig. [3]). Instead we intro- 
duce a measure for the compactness by considering the 
enclosed mass within ellipsoids with the semiaxes a = x, 
h = X and c = x/2. (Note that the hypermassive ob- 
ject is highly deformed. The defined ellipsoids describe 
approximately isodensity surfaces.) Fig. [5] clearly illus- 
trates the different remnant structures of the three mod- 
els based on the LS EoS and confirms our explanation 
given above. The enclosed mass is computed about 8 ms 
after the merging, where the merger time is defined as the 
moment when the gravitational- wave amplitude becomes 
maximal. 
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Basically the same picture as for the LS EoS holds for 
the simulations using the Shen EoS. The luminosity spec- 
tra are displayed in Fig. [6l The generally lower values of 
/peak are explained by the fact that the Shen EoS leads 
to less compact objects than the LS EoS. The differences 
in the peak frequencies are only about 50 Hz (see Tab.|I|. 
Again, assuming Fth = 1.5 results in a higher oscillation 
frequency, while using Fth = 2 yields a lower peak fre- 
quency compared to the model with the fully consistent 
temperature description. The simulations with the Shen 
EoS lead to smaller differences because, as mentioned in 
Sect. El generally lower temperatures are obtained with 
this EoS. Therefore, thermal effects are not as impor- 
tant as in the case of the LS EoS, and an approximate 
inclusion has less influence. 

Our results arc also consistent with the findings sum- 
marized in Sect. |TT1 where we mentioned that the as- 
sumption of perfectly efficient cooling of the neutron star 
merger remnant leads to a more compact remnant struc- 
ture. For this reason, as a consequence of the miss- 
ing thermal pressure support, the peak frequencies were 
found to be higher in Q. 

Note that the low-frequency domain of the 
gravitational-wave luminosity spectra is mainly de- 
termined by the inspiral phase, where thermal effects are 
not important. This explains why the spectra in Figs. 0] 
and [5] are very similar at low frequencies. 



is formed only after 20.8 ms. The difference is a conse- 
quence of the higher compactness of the central object 
in the Fth = 1.5-run. Simulating thermal effects with 
an ideal-gas index of Fth = 2 leads to a delay time of 
21.0 ms, which is in good agreement with the calculation 
employing the full temperature-dependent EoS. We as- 
sume that this consistency is a consequence of the fact 
that the structure of the inner region of the remnant, 
where the collapse sets in, is better reproduced by the 
Fth = 2 choice (see Fig. [S]). The uncertainties of Tdciay 
are important findings because they mean that, depend- 
ing on the specific choice of Fth, the delay time as a very 
basic property of the merger remnant may be incorrect 
by a factor of about 2 in simulations using an ideal-gas 
ansatz. 

Because of the high computational costs of our long- 
term calculations, we cannot determine the delay times 
for our models based on the Shen EoS. The merger rem- 
nants in these models do not collapse to black holes dur- 
ing the simulations. This result is clear, because the Shen 
EoS leads to less compact central objects and it supports 
higher neutron star masses in comparison to the LS EoS. 
Both delays the formation of a black hole. In Tab. |T] we 
give lower bounds on the lifetime of the hypermassive 
remnants for the Shen EoS. 



C. Torus masses 



B. Delay time to black hole collapse 

Besides influencing the oscillation frequency of the 
postmerger remnant, the different remnant structure also 
affects the lifetime of the central hypermassive neutron 
star. Since the remnant rotates differentially, more mass 
can be stabilized than for a rigid rotator. It therefore 
collapses to a black hole after angular momentum redis- 
tribution. Typically, the delay time Tdoiay between the 
merging and the collapse is of the order of several 10 ms 
and depends strongly on the total mass of the system and 
the EoS [il,0,[23- 

The delay time Tdeiay is an important observational 
property, because it determines the maximal length of 
the gravitational- wave emission by the stable postmerger 
remnant. The signal carries information on the EoS of 
high-density matter. The gravitational waves produced 
by the relic black hole and the surrounding torus are of 
fundamentally different character, e.g. at much higher 
frequencies [y, [3 • Finally, the delay time itself may be 
an interesting measure as it depends sensitively on the 
binary system parameters and the EoS. Consequently, an 
accurate determination of this property by simulations is 
required to find out more about the nature of neutron 
star mergers and of supranuclear matter. 

The collapse to a black hole takes place after Tdoiay = 
10.3 ms in the case of our model with Fth = 1-5 and the 
cold LS EoS (see Tab. [l|. Remarkably, in the simula- 
tion using the full temperature dependence a black hole 



Besides gravitational waves, another characteristic fea- 
ture of neutron star mergers of observational relevance is 
the amount of matter forming a torus around the black 
hole emerging from the collapsing central object. These 
systems might be the sources of short gamma-ray bursts, 
provided the energy release from the black holc-torus 
configuration launches a relativistic outflow producing 
gamma-rays by internal shocks. In the context of this 
scenario the torus mass is a crucial quantity, because it 
determines the possible energy release. This torus mass 
depends on the binary parameters and the EoS 0, 0, HO] ■ 

The torus mass can be estimated by the requirement 
that the specific angular momentum of torus matter must 
be larger than the value at the innermost stable circular 
orbit, i.e. the torus mass can be roughly determined 
before the central object collapses into a black hole. A 
description of the iterative procedure is given in Q . The 
formation of the black hole cannot be followed by our 
code, which is why we monitor the amount of matter 
fulfilling the torus criterion as a function of time and 
determine the torus mass shortly before the onset of the 
gravitational collapse. 

The estimated torus masses are listed in Tab. HI The 
values increase towards the end of the simulations, be- 
cause the merger remnants are not yet stationary. In all 
six simulations we observe a continuous contraction of 
the central object. 

Only in the calculations with the LS EoS we can de- 
termine the final torus mass, because only in these cases 
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the gravitational collapse of the merger remnant occurs 
within the simulation time. The torus mass is influenced 
by two effects, which can be disentangled by monitoring 
the temporal evolution of the amount of matter fulfilling 
the torus criterion. After the initial plunge the former 
cores of the stars bounce and reexpand to contract again 
several times until these oscillations are damped after 
some cycles. In the first expansion phase the more com- 
pact remnant of the Fth = 1.5 model leads to a strong 
initial increase of the amount of potential torus matter, 
while this effect is strongly suppressed in the Fth = 2 sim- 
ulation. The fully consistent model with an intermediate 
compactness shows a moderate growth during the first 
decompression stage. (We find the same effect compar- 
ing the T = cases and the consistent models discussed 
in In the subsequent evolution the amount of mat- 
ter with sufficient angular momentum to fulfill the torus 
criterion steadily increases with roughly a constant, but 
lower and similar rate in all models. Therefore, the delay 
time until black hole formation affects sensitively the fi- 
nal torus mass. Although starting with an initially lower 
amount of matter obeying the torus criterion, the model 
with Fth = 2 can therefore end up with a higher final 
value compared to the Fth = 1-5 case, because the latter 
model collapses to a black hole relatively early. The high- 
est final torus mass is obtained for the simulation with 
the full temperature dependence of the EoS, because the 
initial amount of torus matter from the first expansion is 
higher than the one of the Fth = 2 calculation, and the 
delay times are comparable, admitting a long phase of 
angular momentum transport to the outer parts of the 
merger remnant. 

For the models based on the Shcn EoS the amount 
of matter fulfilling the torus criterion 16 ms after the 
merging is given in Tab. HI Note that these values are only 
very crude estimates for the final torus masses, because 
this quantity depends on the delay time, which cannot 
be determined in our simulations with the Shen EoS. 

For both EoSs the evolution of the torus masses in the 
simulations with Fth = 2 follows the models with the 
fully consistent temperature treatment more closely. 

Furthermore, we analyzed the amount of matter that 
becomes gravitationally unbound from the merger site. 
This is of interest for the question how neutron star 
merger events might contribute to the galactic abun- 
dances of heavy elements, e.g. by the production of rapid 
neutron capture nuclei in the ejecta [soj . We find dis- 
agreement in the ejecta masses of about a factor of 2 
between our reference models and the calculations with 
the simplified EoS treatments. 



VII. CONCLUSIONS 

We performed simulations of neutron star mergers 
with microphysical EoSs employing an approximate de- 
scription of thermal effects by an ideal-gas extension of 
the zero-temperature EoSs on the one hand, and us- 



ing the full temperature dependence of the EoSs on the 
other hand. This allows us to assess the limitations of 
the approximate approach and in particular to quan- 
tify the deviations occurring in observational quantities 
like the gravitational-wave signal. For characterizing the 
gravitational-wave emission of the postmerger remnant, 
which is the phase of the event influenced most by ther- 
mal effects, wc determine the dominant gravitational- 
wave frequency of the oscillations of the forming hyper- 
massive object. It is found that the simplified inclusion 
of thermal effects yields values that differ by about 50 to 
250 Hz from the value obtained in the simulation with 
the fully consistent temperature treatment. The differ- 
ences can be understood as a consequence of the different 
remnant compactness, which is caused by the different 
strengths of the thermal pressure contributions. 

Furthermore, the simplified approach affects sensi- 
tively the delay time between merging and the formation 
of a black hole. Again, the discrepancies can be explained 
by the different remnant structures, where a more com- 
pact object collapses to a black hole earlier. The lifetime 
of the hypcrmassive remnant can deviate from its true 
value by a factor of 2 in simulations that use an ideal-gas 
like ansatz to describe thermal effects. These uncertain- 
ties are in particular important, because they influence 
also the duration of the gravitational- wave emission from 
the postmerger remnant and thus the integrated power 
in the relevant frequency bands, which provide informa- 
tion about the high-density EoS. In addition, errors in 
the delay time also have an effect on the estimates of 
the amount of matter remaining outside the black hole 
and the forming accretion torus after the gravitational 
collapse of the central object. 

Our results show that estimates of this disk mass by 
means of simulations that use the ideal-gas approxima- 
tion, are uncertain by up to 30% of the full temperature 
result. These uncertainties arise because of the effects of 
the approximate EoS treatment on the delay time as well 
as the remnant structure and evolution. 

The comparison of simulations with the LS EoS sug- 
gests that Fth = 2 might be a better choice in particu- 
lar to model the delay time, although an agreement in 
all relevant quantities cannot be achieved. For instance, 
considering the frequency of the postmerger ringdown, 
Fth = 1.75 may yield more accurate results. However, 
conclusions based on such a small number of models (with 
respect to the investigated EoSs and the binary configu- 
rations) have only a limited significance, for which rea- 
son we refrain from specifying optimal values for different 
needs here. 

In summary, we stress that the interpretation of re- 
sults from simulations using the approximate implemen- 
tation of thermal pressure contributions should take into 
account the uncertainties presented in this study (see 
Tab. IJ). This has an impact in particular on attempts 
to constrain the high-density EoS by observational data 
from merging events. 

In addition, we expect that for more massive binary 
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configurations the deviations due to the simplified treat- 
ment even increase, because mergers of such systems 
yield higher temperatures [1, 01 • 

We note that our findings apply also to other than 
the investigated zero-temperature EoSs when an ideal- 
gas ansatz is used to account for temperature effects. For 
instance, it was shown that picccwisc polytropic EoSs can 
be chosen in order to match microphysical zero temper- 
ature EoSs [3l|. Also in this case the thermal behav- 
ior of the EoS needs to be implemented in an approxi- 
mate manner giving rise to the associated uncertainties. 
Therefore, we conclude that microphysical EoSs that in- 
clude thermal effects consistently are highly demanded 
for astrophysical applications as the ones discussed here, 
in particular if a high quantitative accuracy is required. 
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